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ABSTRACT 

We study the properties of X-ray galaxy clusters in four cold dark matter models with 
different baryon fraction CIbm ranging from 5 to 20 per cent. By using an original 
three-dimensional hydrodynamic code based on the piecewise parabolic method, we 
run simulations on a box with size 64 h~ l Mpc and we identify the clusters by se- 
lecting the peaks in the X-ray luminosity field. We analyse these mock catalogues by 
computing the mass function, the luminosity function, the temperature distribution 
and the luminosity-temperature relation. By comparing the predictions of the different 
models to a series of recent observational results, we find that only the models with 
low baryonic content agree with the data, while models with larger baryon fraction are 
well outside the l-cr errorbars. In particular, the analysis of the luminosity functions, 
both bolometric and in the energy band [0.5-2] keV, requires £Ibm~ 0.05 when we fix 
the values h = 0.5 and n = 0.8 for the Hubble parameter and the primordial spectral 
index, respectively. Moreover we find that, independently of the cosmological scenario, 
all the considered quantities have a very little redshift evolution, particularly between 
z = 0.5 and 2 = 0. 

Key words: Hydrodynamics - Large-scale structure of the Universe - X-ray: galaxies, 
general - Dark matter - Galaxies: clusters -Cosmology: theory 
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1 INTRODUCTION 

Galaxy clusters are the most extended gravitationally bound 
systems of the Universe. In spite of their size, they are simple 
objects: their dynamics is governed essentially by the bal- 
ance between gravitational and pressure forces, controlled by 
shock heating and adiabatic compression. Other processes, 
like radiative cooling in cluster cores or population III stars 
heating, may play a significant role only in some aspects of 
their history. For these reasons, the evolution and present 
day properties of clusters are very sensitive to the funda- 
mental cosmological parameters and to the initial power 
spectrum. Therefore they provide an ideal tool to study 
the formation of the structures of the universe. Analytic 
techniques, as the Press- Schechter (1974) formalism, and 
numerical simulations, based on the Zel'dovich (1970) ap- 
proximation or on more accurate N-body codes, have been 
extensively used to put strong constraints on the cosmo- 
logical scenarios by comparing the model predictions to the 
observed abundances and clustering properties (see e.g. Bor- 
gani et al. 1997 and references therein). 

It is well known that clusters have a large X-ray emis- 
sion, essentially produced by thermal bremsstrahlung. The 



quality of the observations in this band have been recently 
largely improved, providing new data which potentially can 
discriminate between different cosmological models. In order 
to compare these data to the numerical simulations we must 
include a reliable treatment of the hydrodynamic quantities, 
such as temperature, pressure and energy. For these reasons 
in the last years there has been a large effort of the cosmo- 
logical community to develop new numerical codes which 
combine the standard N-body techniques, which solve the 
dynamics of the dark matter component, to hydrodynamic 
methods for properly following the evolution of the bary- 
onic component (Hernquist & Katz 1989; Evrard 1990; Cen 
1992; Steinmetz & Muller 1993; Ryu et al. 1993; Bryan et 
al. 1995; Gnedin 1995; Pen 1997; for a comparison of the dif- 
ferent algorithms see Kang et al. 1994b). These codes have 
been applied for simulating in various cosmological models 
galaxy cluster properties and evolution and for calculating 
their masses, X-ray emission, mean temperature and the re- 
lations between these quantities, that is the fundamental in- 
formation that is usually obtained from observational data 
(see e.g. Kang et al. 1994a; Bryan et al. 1994a,b; Navarro, 
Frenk & White 1995; Pen 1996; Cen 1997). 
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In this paper we apply a hydro dynamic code based on 
the coupling of the Piecewise Parabolic Method with the 
Particle-Mesh N-body code (Gheller, Pantano & Moscar- 
dini 1997) to study the large-scale distribution of the X-ray 
clusters in the framework of critical-density cold dark mat- 
ter (CDM) models with a high baryon content (up to 20 
per cent). This kind of model, allowed by the present un- 
certainties in the determination of the baryon density Q,bm 
of the universe, when considered with a small tilt of the 
primordial spectral index (n ~ 0.8), has been suggested as 
possible solution of the problems of the standard (i.e. with 
flflM ~ 0.05 and n = 1) CDM model, namely the high 
small-scale power and the low first peak in the microwave 
anisotropy power spectrum (White et al. 1996). For a dis- 
cussion about tilted models, see also Lucchin & Matarrese 
(1985); Vittorio, Matarrese & Lucchin (1988); Liddle & Lyth 
(1993). 

The plan of the paper is as follows. In Section 2 we in- 
troduce the cosmological models that we consider. In Section 
3 we present the numerical code and the method of cluster 
identification. The results are shown in Section 4. The main 
conclusions are drawn in Section 5. 



2 COLD DARK MATTER MODELS WITH 
HIGH BARYON CONTENT 

For a long period the CDM scenario has been the reference 
model for the interpretation of the observational data on the 
large scale of the universe. Its standard version assumes a 
flat universe with a density parameter Slo = 1, a Hubble 
constant h = 0.5 (in units of 100 km s _1 Mpc -1 ), a baryon 
contribution to the density Q.bm fixed by the standard the- 
ory of big bang nucleosynthesis (BBN), and primordial fluc- 
tuations with Gaussian distribution and power spectrum 
P(k) oc k" , with n = 1. However, the normalization im- 
plied by the COBE detection of the microwave anisotropics 
(Smoot et al. 1992; see also Bennett et al. 1996) changed the 
situation, giving a spectrum with too much power on scales 
smaller than 10 hT 1 Mpc. As a consequence, the CDM model 
is not able to reproduce either the clustering properties of 
galaxies and the distribution and abundances of clusters. 

In order to solve the CDM problems, many alternatives 
have been proposed and discussed to reduce the short-scale 
power: tilted CDM models, i.e with a primordial spectral 
index n < 1; mixed dark matter models, i.e. with the addi- 
tion of a hot component corresponding to approximately 20 
per cent of the total density; open CDM models, i.e. with 
fio < 1; Lambda CDM models, i.e. with a non- vanishing 
cosmological constant fixed such as its contribution to the 
density is Qa = 1 — fio; CDM models with a small Hubble 
constant, i.e. with h < 0.5. For a general review on the CDM 
model and its variants, see Liddle & Lyth (1993) and Coles 
(1996). 

Very recently, it has been suggested (White et al. 1996) 
that also an increase of the baryon content £Ibm can pro- 
duce critical-density CDM models which fit the observa- 
tional data in a reasonable way. In the past, the mean baryon 
density has been considered as a fixed quantity for most of 
cosmological models, in agreement with the predictions of 
the standard theory of BBN and the local estimates of light 



elements abundances: Q.BMh 2 = 0.0125 ± 0.0025, where the 
uncertainty is at 95 per cent confidence level. 

The possibility of varying Q.bm has been introduced 
because some recent observations have suggested that its 
estimate is much more uncertain than previously thought. 
In fact the first measurements of the deuterium abundance 
at high redshift yield very discrepant values. Small Osm, 
in good agreement with the BBN predictions, are implied 
from different measurements of the neutral hydrogen column 
density in high-redshift clouds (Rugers & Hogan 1996a,b; 
Songaila, Wampler & Cowie 1997). On the contrary, the 
analyses performed by Tytler, Fan & Buries (1996) suggest 
Q.BMh 2 = 0.024 ± 0.002 ± 0.002 ± 0.001, with the la uncer- 
tainties being statistical, systematic and theoretical, respec- 
tively. These last results are also consistent with the lower 
limits obtained by Rauch et al. (1997) by comparing the ob- 
served flux decrement distribution function from a sample 
of seven high resolution QSO spectra to simulations of the 
Lya forest. Similar conclusions have been reached also by 
Weinberg et al. (1997) by computing the amount of neutral 
hydrogen present in the high-redshift intergalactic medium 
necessary to produce the Lya absorption in the QSO spec- 
tra. Finally, the determination of the baryonic fraction in 
clusters of galaxies (White et al. 1993; White & Fabian 
1995; Elbaz, Arnaud & Bohringer 1995) leads, in a critical- 
density universe, to high values of Q,bm- For example White 
& Fabian (1995) have calculated: 

Q bm _ nl4 +o.Q8 ( h y^ 2 

again at the 95 per cent confidence. 

By using a semi-analytical approach, White et al. (1996) 
carried out a general exploration of CDM models with high 
fisM, allowing also the Hubble parameter h and the spec- 
tral index n to vary. The predictions of these scenarios 
have been compared with the observations of the clustering 
properties of galaxies, the cluster abundances, the statis- 
tics of the peculiar velocities, the formation of high-redshift 
objects (i.e. damped Lya systems, Lyman break galaxies, 
quasars and clusters) and the cosmic microwave background 
anisotropics. Their final suggestion is that the models with 
f2sM in the range [0.1-0.2], a Hubble parameter h w 0.5 and 
a small tilt in the primordial spectrum (n « 0.8), are in good 
agreement with all these data. A high baryonic fraction in 
fact helps "naturally" to suppress the short-scale power, as 
baryon collapse is stopped till decoupling, and, at the same 
time, amplifies the first peak of the CMB spectrum, compen- 
sating for the loss of height introduced by the tilt (needed 
to avoid too low values of h, excluded by the observational 
data). At the end of their analysis they conclude that these 
models can represent at the moment a viable alternative in 
the framework of critical-density CDM models. 

In this paper we consider models where the universe is 
spatially flat, i.e. Qo — Qdm + Qbm = 1. Here Qdm is 
the density parameter of the (cold) dark matter component. 
We assume a vanishing cosmological constant (Qa = 0). 
The initial spectrum of perturbations is defined as P(k) oc 
k n T 2 (k), where T(k) is the CDM transfer function. We use 
the expression obtained by Bardeen et al. (1986): 

TM . Mi_^) x , 1+ , 89 , + 
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Table 1. The parameters of the cosmological models. Column 2: 
the density parameter Ho; Column 3: the baryon density Qbm\ 
Column 4: the primordial spectral index n; Column 5: the Hubble 
parameter h; Column 6: the shape parameter T; Column 7: the 
spectrum normalization erg. 



Model 


n 


Qbm 


n 


h 


r 


o"S 


BM05 


1.0 


0.05 


0.8 


0.5 


0.44 


0.77 


BM10 


1.0 


0.10 


0.8 


0.5 


0.41 


0.72 


BM15 


1.0 


0.15 


0.8 


0.5 


0.39 


0.66 


BM20 


1.0 


0.20 


0.8 


0.5 


0.36 


0.61 



(16.1 9 ) 2 + (5.46g) 3 + (6.71q) 4 ]- 1/4 , (2) 

where q = k/hY. The shape parameter F takes into account 
the dependence on the Hubble parameter h, on the total den- 
sity Oo and on the baryon density Qbm (Sugiyama 1995): 

T = n /iexp(-a SM - Qbm/^o) ■ (3) 

As suggested by White et al. (1996), we have fixed 
h = 0.5 and n = 0.8, allowing the baryon content to vary. 
We consider four different values for Obm: the usual value 
determined by the standard theory of BBN (Qbm = 0.05); 
a baryon abundance consistent with the low deuterium 
measurements and with the more recent BBN calculations 
(Qbm = 0.1); a value close to the estimate resulting from 
the cluster baryon fraction in the case of critical-density uni- 
verse {O.BM = 0.15); a more extreme case representing the 
upper limit of the range of observed cluster baryon fraction 
(fisM = 0.2). In the following we will label these four models 
as BM05, BM10, BM15 and BM20, respectively. The nor- 
malization of the spectrum, usually parametrized by as, the 
matter rms fluctuation in a top-hat sphere of radius 8 h^ 1 
Mpc, is defined by the four-year COBE data (Bunn & White 
1997). The cosmological parameters used for the different 
models are summarized in Table 1. Since the power at small 
scales (i.e. on large k) increases when Qbm decreases, we 
can expect a faster evolution in the models with low baryon 
content, with a larger production of big mass overdensities. 

3 NUMERICAL METHOD 
3.1 The code 

The simulations have been evolved by using the numeri- 
cal code presented in Gheller et al. (1997; see also Gheller, 
Moscardini & Pantano 1996), where we give a complete de- 
scription of the method and we show the results of various 
numerical tests. Here we only summarize the main charac- 
teristics of the code. 

The hydrodynarnical part has been developed by us- 
ing the Eulerian version of the Piecewise Parabolic Method 
(PPM; Colella & Woodward 1984) which ensures at least 
second-order (up to the fourth-order, in the case of smooth 
flows and small timesteps) accuracy in space and second- 
order accuracy in time. The high accuracy of this method 
allows minimization of errors due to the finite size of the 
cells of the grid and leads to a spatial resolution close to the 
nominal one (i.e. one grid). In a cosmological framework, 



the basic PPM technique has been modified to include the 
gravitational interaction and the expansion of the universe. 
Particular care has been devoted to the calculation of the 
gas internal energy and a double formulation of the energy 
equation has been used in order to avoid large errors in 
the computation of thermodynamical quantities when the 
kinetic energy is very large compared to the internal one. 

The hydrodynarnical part has been coupled to a Particle 
Mesh (PM) N-body code (Hockney & Eastwood 1981) that 
describes the evolution of the dark component. The standard 
PM code has been modified in order to allow non-constant 
timesteps equal to those used in the integration of the hydro- 
dynamical equations. This is obtained by replacing the stan- 
dard second-order leapfrog method by a second-order two- 
step Lax-Wendroff scheme. Densities and forces are com- 
puted by using the cloud-in-cell interpolation scheme. The 
coupling is obtained by calculating, by the usual FFT pro- 
cedure, the gravitational field due to both components. 

In the simulations presented in this paper we neglected 
the atomic processes for radiative cooling since the cooling 
time of the hot gas produced in clusters is longer than a Hub- 
ble time. Moreover the only considered heating processes 
are the adiabatic compression and the entropy generation 
at shock fronts. 

For each cosmological model we run one simulation with 
the initial conditions given by the same random sequence. 
The initial redshift, fixed in such a way that the maximum 
initial density fluctuation is less than unity, is approximately 
z ~ 20 for all models. The box-size has been fixed to 64 h^ 1 
Mpc and the number of computational cells is 128 3 . Conse- 
quently the nominal spatial resolution, which our numerical 
tests have shown to be very close to the effective one (Gheller 
et al. 1997), is 0.5 h~ Mpc. In the analysis of the results we 
have to be careful to properly evaluate the spurious effects 
that the limited box-size and the finite grid resolution can 
produce on the results. In fact the grid resolution prevents 
us from following the behaviour of matter inside a cell el- 
ement. This can lead to underestimation of quantities like 
density and temperature, with a direct influence on the evo- 
lutionary history of the X-ray clusters. On the other hand, 
the limited size of the box has the effect of suppressing the 
large-scale power, reducing the possibility of forming very 
bright X-ray clusters with luminosity larger than 10 45 erg 
s _1 . This kind of object is also likely to be partially missed 
in our simulations because it is rare: observational data (e.g. 
Henry & Arnaud 1991; Ebeling et al. 1997) show that in a 
box of 64 h^ 1 Mpc one expects to find at most one such 
object. 



3.2 Cluster Identification 

The first step of the data analysis consists of the identifi- 
cation of the X-ray clusters. We first calculate the emissiv- 
ity due to thermal bremsstrahlung in a fully ionized plasma 
(X = 0.76, Y = 0.24) with temperature T (Rybicki & Light- 
man 1979; see also Evrard 1990): 

e„ = 6.83 x 10~ 3S Z 2 n e mT' 1/2 x 

g(hu/kT) e- hv/kT erg s" 1 cnT 3 Hz" 1 , (4) 
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where g(hv/kT)e~ hL/ < kT is an average Gaunt factor, Z is 
the charge number, m and n e are the ion and the electron 
density, respectively. 

The bolometric emissivity, using a unit Gaunt factor, is 
defined as: 

tbol = / dve v 

Jt) 

= 1.42 x KT 27 Z 2 n;n e T 1/2 erg s~ x cm~ 3 . (5) 

Then the energy radiated within a given energy band E\ - 
E2 can be expressed as 

tband = fband(T) €bol , (6) 

where 

rE 2 /kT 

f band (T) = dr, g(n) e"" . (7) 

J Ei/kT 

The band limited X-ray emission L x from a given volume is 
computed integrating the previous expression over the rele- 
vant volume. Using the discretization of the simulation, the 
X-ray luminosity is 

Lx = 1.25 x 10- 27 m- 2 ^2pl M T^ 2 f band (Ti) erg s" 1 , (8) 

where the sum runs over cells within the volume, p BM is 
the baryon density and m p is the proton mass. Since we 
are interested only in X-ray emitting regions, it is safe to 
assume that the gas is completely ionized. We also set the 
Gaunt factor g = 1.2: this gives an accuracy of ~ 20 per 
cent for the results. 

At this point, in order to identify the clusters, we select 
the cells with L x > 10 40 erg s _1 which are also local maxima 
in the X-luminosity field (i.e. their X-luminosity is greater 
than that of the 26 neighboring cells). These identify the 
cluster centres. A cluster is defined as the sum of the centre 
plus the 26 surrounding cells. In this way the total volume of 
a cluster equals the volume of a sphere of comoving radius 
0.93 hT 1 Mpc, as appropriate for present observed X-ray 
clusters. In order to avoid double counting of the cells, the 
distance between cores is checked: if two cores are closer than 
2 h^ 1 Mpc the fainter cluster is rejected from the catalogue. 
Finally, the cluster luminosity and the mass are calculated 
as the sum of the luminosity and density of each of its cells, 
respectively, while the temperature is defined as the average 
over the whole cluster volume. 



4 RESULTS 

4.1 Global properties 

In Figure 1, we show a snapshot of the results in a slice 
of 64 x 64 x 0.5 h~ 3 Mpc 3 at z = for each of the four 
models. The baryonic matter density field Q BM , the dark 
matter density field Q DM , the gas temperature T and the X- 
ray emission L x are presented in Figures la, lb, lc and Id, 
respectively. Since the initial spectrum of fluctuations of all 
simulations has been generated by using the same random 
number sequence, the positions of the final structures are 
quite similar and the densities of each component scale ap- 
proximatively according to their mean cosmic values. Matter 



concentrates on filamentary structures and clusters form at 
the intersection of several filaments. 

In the \ow-Qbm models, because of the larger power on 
small scales in the initial spectrum and the lower background 
pressure, shocks form earlier and are stronger than in the 
high-^BM models. In the latter case, on the other hand, 
we can observe a stronger X-ray emission due to the higher 
baryon content of these models. 

A close-up at four different redshifts of the most lu- 
minous X-ray cluster found in the BM05 simulation is pre- 
sented in Figure 2. The slice is 20 x 20 x 0.5 /i~ 3 Mpc 3 . It is 
evident the ongoing process of merging and the final virial- 
ized state characterized by an extended central isothermal 
region at high temperature. The cluster tends to a spherical 
geometry and deviations from sphericity are due to the clus- 
ter memory of its merging history (see e.g. Tormen 1997). 
These satellites, however, being faint X-ray sources, may not 
be observed in a X-ray map. Similar comments can be re- 
peated for the other models (not shown here). 

Dark matter structures appear typically less concen- 
trate than the baryonic counterpart both in clusters and 
in filaments. In Table 2 we present the rms of the DM and 
BM density fields (p DM and o BM respectively) computed 
on the cell-size scale (0.5 7i~ Mpc) and the mean tempera- 
ture < T > (in Kelvin degrees) at redshifts z = 1 and z = 0. 
The density contrasts of the two components are normalized 
to the corresponding mean cosmological values. A detailed 
analysis of the time evolution of the rms shows that at very 
high redshifts, due to absence of pressure forces, DM col- 
lapses faster than the baryonic counterpart. However, in all 
of our models, starting from about z ~ 3, baryons tend to 
concentrate more than the dark component. The absence of 
dissipative phenomena produces a spreading of DM around 
the minima of the gravitational potential, while BM, which 
tends to thermalize, concentrates there. As time goes by, 
also dark particles fall toward the centre of the potential 
well and at the final time the rms of the two components is 
quite similar (see Table 2). 

Comparing the different models, we observe that both 
a BM and o DM decrease with the increasing of the baryon 
fraction both at z = 1 and z — 0. This behaviour is mainly 
related to the characteristics of the initial spectrum in the 
different models. Furthermore structures virialize earlier in 
low-f2sM models and shocks are stronger and the final tem- 
peratures higher than in the high-QsM models, because of 
the lower pressure and higher densities present initially in 
these models. Between 2 = 1 and z — there is a slightly 
faster evolution of structures in the case of a high baryon 
content and this is probably favoured by the lower temper- 
atures produced in this case. 

4.2 X-ray cluster mass 

Clusters in our simulations are identified through their X- 
ray emission and their characteristic properties, like the total 
mass, the total luminosity and the mean temperature, are 
computed by integrating or averaging over a fixed number of 
cells, as explained in Section 3.2. The values of these quan- 
tities could then be affected by the choice of the number of 
cells used for the calculation. Several tests have shown that 
our procedure provides good estimates for the temperature 
and the luminosity of the cluster. In fact the temperature is 
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Table 2. The global properties of the model simulations. The rms 
of the dark matter, cr DM , the rms of the baryonic matter, c BM , 
and the mean temperature < T > (in Kelvin degrees), computed 
on the cell scale, are shown at redshift 2 = 1 (Columns 2, 3 and 
4) and at z = (Columns 5, 6 and 7). 



Table 3. The X-ray emissivity (in units of 10 40 ergs s _1 ft 3 
Mpc -3 ) for the gas (jgas) and for the clusters (j c i) at various 
redshifts for the different models. 



Models 


a DM 




< T > K 


a DM 




< T > K 






2 = 


1 




2 = 





BM05 


3.27 


3.40 


4.86 x 10 5 


7.31 


7.34 


1.04 x 10 6 


BM10 


2.80 


3.14 


3.86 x 10 5 


6.97 


7.00 


9.24 x 10 5 


BM15 


2.44 


2.73 


2.97 x 10 5 


6.05 


6.27 


7.68 x 10 5 


BM20 


1.94 


2.20 


2.20 x 10 5 


5.33 


5.45 


6.35 x 10 5 



BM05 BM10 BM15 BM20 



Jgas Jci Jgas J cl Jgas J cl Jgas Jcl 

2 = 1 0.17 0.12 0.49 0.39 0.66 0.49 0.73 0.55 

2 = 0.7 0.21 0.16 0.54 0.43 0.79 0.62 0.98 0.79 

2 = 0.5 0.20 0.15 0.63 0.46 0.88 0.73 1.53 0.86 

2 = 0.2 0.19 0.14 0.51 0.39 0.87 0.73 1.29 1.12 

2 = 0.16 0.13 0.59 0.50 0.96 0.81 1.05 0.85 



almost uniform over regions greater than those over which 
we do our averaging. Moreover the luminosity is propor- 
tional to the square of the baryonic density and then its 
value depends essentially on the cells with highest density 
which represents the centre of the cluster and which are con- 
sequently always included in our integration. The inclusion 
or missing of some low-density cells does not affect sensi- 
bly our estimates for the temperature and the luminosity, 
but could affect the calculation of the cluster total mass. In 
summing over a fixed number of cells we tend to overesti- 
mate the mass of objects which are not as extended as our 
reference volume (27 cells) and which are also usually char- 
acterized by density lower that the mean cluster values. The 
opposite is true for large clusters, although in this case we 
have verified that, for our choice of the reference volume, the 
error introduced in the estimation of the mass is less severe 
and it is at most a factor two. Therefore we have restricted 
our analysis to clusters with mass greater than 10 14 M Q (for 
h = 0.5). This lower limit has been fixed in order to avoid the 
inclusion of too low-mass objects which can be affected by 
a large error in the mass estimate and whose properties dif- 
fer considerably from those of a typical galaxy cluster (e.g., 
their mean density is much lower than that determined from 
observations). Furthermore their spatial distribution could 
be affected by the X-luminosity selection criterion used to 
build our cluster catalogue. 

In Figure 3 we present, for the four models, the num- 
ber Nm of clusters with mass greater than 10 14 Mq found in 
the whole simulation box at various redshifts. The number 
of clusters decreases with increasing Qbm at any redshift. 
This is a consequence of the amount of power on small scales, 
which decreases as we increase the baryonic fraction. In the 
high-f^BM cases the collapse of massive objects occurs later 
and, then, if we select clusters by their mass, the models 
with high baryonic fraction presents less clusters than low- 
Osm models. The behaviour of the mass function depends 
on the distribution of the total (dark plus baryonic) density 
and then it is not the preferred quantity for discriminating 
between our models. Other quantities, like the X-ray lumi- 
nosity, that depends directly on the BM density, appear to 
be a more useful quantity for the problem studied in the 
present paper. 

Notice that the estimated Nm roughly agrees with the 
observational data for all the models. This is an expected 
result, because of the spectra normalizations. In fact it is 
known that the present abundance of galaxy clusters re- 



quires, in the framework of the critical-density models and 
almost independently of the shape of the primordial spec- 
trum, a normalization a$ m 0.6, with a quite large uncer- 
tainty (e.g. White, Efstathiou & Frenk 1993; Viana & Liddle 
1996; Eke, Cole & Frenk 1996). 

4.3 X-ray cluster luminosity 

In Figure 4 we present the number of clusters with luminos- 
ity greater than L x = 10 42 and 10 43 erg s _1 (A/42 and -/V43, 
respectively). This luminosity-based selection criterion leads 
to results opposite to those obtained by selecting clusters by 
their masses. In fact at any epoch the number of clusters in- 
creases with increasing Qbm- Furthermore the evolution of 
N42 and A^43 with redshift is not monotonic: the number 
of clusters tends to grow until a turn-around redshift after 
which it starts to decrease. This behaviour is common to all 
the models and for both the adopted minimum luminosity. 
The only exception is ^43 in the case of BM05, which grows 
continuously with time, but, because of the small number, 
this might be not statistically significant. The turn-around 
is due to the balance between the mechanisms driving the 
cluster evolution, and it is an indication of the epoch when 
the merging processes of different structures start to dom- 
inate over the gravitational collapse of each single object. 
In fact the merging leads to larger but smoother structures. 
Since the X-ray emission is proportional to the square of the 
baryonic density, lower luminosities are expected. The effect 
of the merging processes is also shown by the simultaneous 
decrease in the number of fainter clusters (not reported in 
the figure). The turn-around redshift becomes lower with in- 
creasing the baryon fraction ranging from 2 = 0.7 for BM05 
to z = 0.2 - 0.5 for BM20. This is due to the delayed evolu- 
tion of the structures in high-fisM models. 

This behaviour is confirmed by the X-ray emissivity per 
unit comoving volume due to both the gas in its entirety, 
jgas, and the clusters, j c i (see Table 3). For each model the 
two quantities evolve in a parallel way, indicating that clus- 
ters emit roughly a constant large fraction of the total X-ray 
radiation. Similar results have been obtained, for the stan- 
dard CDM model, also by Kang et al. (1994a) and Bryan et 
al. (1994a). 

The fundamental differences between the cluster abun- 
dances in mass (Nm) and in X-ray luminosity (N42 and A/43) 
must be kept in mind when the results are compared to the 
observations. In Figure 5 we present the luminosity function 
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for the four models computed at five different redshifts, inte- 
grated over the whole range of frequencies. Very bright clus- 
ters with luminosity greater than 10 45 erg s _1 are missing in 
our simulations. The lack of such clusters can be related to 
two different effects. Firstly, the size of our computational 
box limits the amount of large-scale power which we can 
follow in the simulation, and consequently the maximum 
temperature that can be produced. Secondly, the grid res- 
olution is likely to underestimate the highest density peaks 
where brightest clusters are expected to form. 

The luminosity functions have been fitted by using a 
two-parameter function: 

n{L)dL = n L' a dL , (9) 

where n(L)dL is the comoving density of clusters with lumi- 
nosity between L and L + dL, no is in units of 10" 6 /i 3 Mpc -3 
and L is in units of 10 44 erg s _1 . The results of the fits 
are presented in Table 4. Notice that for our data a two- 
parameter fit is more appropriate than the usual Schechter 
function, as in our results the expected bend at high L is 
not present, for the reasons given above. 

The values of the slope a cannot discriminate between 
the four models, as the differences are in general within the 
1-a errorbars. The parameter a shows for all the models 
a slight negative evolution with redshift from z = 0.5 to 
z = 0; this behaviour is thought to be the effect of the 
ongoing processes of gravitational collapse, that produces 
more concentrated and bright structures, and merging, that 
leads to the formation of larger and more massive objects 
at the expenses of the smaller ones. For the normalizations 
no, which is the comoving number density of objects with 
luminosity equal to 10 44 erg s _1 , we can make considerations 
similar to those previously made for N42 and N43. 

The hatched region in the figure shows the observational 
data (with 1-a errorbars) of Ebeling et al. (f997) which re- 
fer to the ROSAT Brightest Cluster sample containing f99 
objects with redshift z < 0.3. The observational curves have 
been fitted by the authors by using a three-parameter func- 
tion: 

n(L) = Aexp(-L/L*)L" Q , (10) 

where A is in units of 10~ 7 Mpc~ 3 (10 44 erg s" 1 )"' 1 and L* is 
in units of 10 44 erg s _1 . The values of the fitting parameters 
are A = 6.41±°;™, L* = 37.2+J 6 8 4 and a = 1.841°;°°. 

We can compare these data with the results of the simu- 
lations at low redshifts. The models with high baryonic con- 
tent (BM15 and BM20) have a luminosity function which is 
significatively too high with respect to the observations. On 
the contrary the BM05 model and (much more marginally) 
BM10 are in better agreement with the data. 

Similar conclusion can be obtained if we consider the 
luminosities in the energy band [0.5-2] keV. The fitting pa- 
rameters for our simulations are reported in Table 5, while 
the comparison with two different observational datasets is 
shown in Figure 6. The vertically hatched region refers again 
to the Ebeling et al. (1997) sample, whose luminosity func- 
tion has been fitted by a three-parameter relation (FloT) with 
A = 3.32±g;f!j, L* = 5.70toja and Q = l-8 5 ±o'.09- The hor- 
izontally hatched region shows instead the results obtained 
by De Grandi (1996) using a complete flux- limited ROSAT 
sample selected from the ESOKP redshift survey: in this 
case the fitting parameters are A = 4.51, L* = 2.63^Q 5g and 



a = 1.32_q 23 . These two determinations of the luminosity 
functions are in good agreement for luminosities larger than 
« 2x 10 43 erg s^" 1 , while for smaller L x the De Grandi (1996) 
results are approximately a factor 3 smaller than the Ebel- 
ing et al. (1997) ones, increasing the discrepancies between 
the observations and the model predictions of the models 
with high Qbm- 

In this energy band there is a further luminosity func- 
tion determined by Burns et al. (1996), always by using im- 
ages from the ROSAT all-sky survey. Because of the large 
errorbars, it completely overlaps both the previous results 
and for clarity we prefer do not show it in Figure 6. How- 
ever, since in this dataset also nearby poor clusters have 
been considered, this result allows to extend the previous 
considerations also to smaller X-ray luminosities (less than 
10 42 erg s _1 ), not included in the other datasets. 

4.4 X— ray cluster temperatures 

In Figure 7 we show the redshift evolution of the distribution 
of the cluster mean temperature for the four models. The 
temperatures have been calculated as emission-weighted av- 
erages because this is the quantity which is also usually es- 
timated from the observations. The absence of clusters with 
temperatures above 4 keV is mainly related to the limited 
size of the box and their rareness. These model predictions 
can be compared with the observations. In the figure the 
hatched region refers to the temperature distribution ob- 
tained by Henry & Arnaud (1991) from a set of local (z « 0) 
clusters: 

n(T) = (l.g±° - 8 B x 10 -3 h 3 Mpc" 3 keV -1 ) T" 4 ' 7 * ' 5 , (11) 

where the uncertainties are 1-a errorbars and T is expressed 
in keV. 

All the models are in quite good agreement with ob- 
servations in the overlapping range. Temperature is in fact 
less sensitive than luminosity to the details of the density 
distribution and it is related to the maximum wavelength 
A of non-linear waves. In fact the post-shock temperature 
is of the order of T oc (if A) 2 , where H is the Hubble con- 
stant. Low-f2_gjvf models have a higher normalization of the 
primordial spectrum and longer wavelengths can reach the 
non-linear regime at the final time producing higher val- 
ues of the temperature. This phenomenon is likely to be 
strengthened by the different mean background pressure of 
the various realizations, which is lower in 1ow-S7sa/ models. 
Both effects could explain the higher number of objects with 
temperature larger than about 1 keV found with decreasing 

Q.BM- 

We have also analysed the redshift evolution of the tem- 
perature distributions in our simulations. We found that be- 
tween z = 1 and z — such distributions are almost con- 
stant for BM05 and BM10 models, with a slow increase in 
the number of high-temperature objects. In these models, by 
z ~ 1 the regions heated up by shocks at almost a uniform 
temperature are larger than the integration volume of our 
cluster identification method. In general, from that moment, 
the temperature of these regions increases because of adia- 
batic compression and merging processes and this explains 
the rise in the number of high-temperature clusters towards 
z = 0. The models BM15 and BM20, instead, evolve rapidly 
between = 1 and z = 0.5, showing in particular a strong 
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Table 4. The parameters of the fits of the X-ray cluster bolometric luminosity function n{L) = tiqL a at various redshifts for the 
different models. 







BM05 


BM10 


BM15 


BM20 






a 


n 


a 


n 


a 


n 


a 


n 


2 


= 1 


1.82± 0.09 


6.73± 0.06 


1.72± 0.08 


18.95± 0.23 


1.59± 0.07 


31.66± 0.38 


1.74± 0.06 


18.56± 0.15 


z 


= 0.7 


1.73± 0.15 


9.72± 0.18 


1.72± 0.09 


21.55± 0.30 


1.71± 0.05 


26.60± 0.24 


1.73± 0.07 


28.44± 0.33 


z 


= 0.5 


1.74± 0.20 


8.42± 0.20 


1.75± 0.16 


19.46± 0.47 


1.76± 0.09 


25.66± 0.36 


1.69± 0.08 


33.20± 0.46 


z 


= 0.2 


1.61± 0.07 


12.50± 0.12 


1.68± 0.04 


20.97± 0.12 


1.69± 0.05 


29.12± 0.24 


1.69± 0.11 


35.59± 0.72 


z 


= 


1.41± 0.11 


23.11± 0.48 


1.67± 0.11 


18.68± 0.31 


1.62± 0.05 


31.69± 0.25 


1.66± 0.03 


33.95± 0.15 



Table 5. The parameters of the fits of the X-ray cluster luminosity function n(L) = ngL a computed in the [0.5-2] keV band at various 
redshifts for the different models. 







BM05 


BM10 


BM15 


BM20 






a 


no 


a 


no 


a 


no 


a 


no 


2 


= 1 


1.47± 0.11 


11.73± 0.17 


1.43± 0.07 


19.26± 0.24 


1.39± 0.10 


19.19± 0.31 


1.40± 0.11 


10.17± 0.13 


2 


= 0.7 


1.58± 0.40 


8.40± 0.49 


1.55± 0.11 


15.73± 0.26 


1.58± 0.11 


12.90± 0.20 


1.45± 0.12 


18.25± 0.32 


2 


= 0.5 


1.81± 0.16 


5.52± 0.04 


1.58± 0.12 


15.09± 0.28 


1.52± 0.07 


22.51± 0.30 


1.57± 0.10 


28.56± 0.48 


2 


= 0.2 


1.44± 0.11 


13.71± 0.27 


1.55± 0.06 


15.87± 0.17 


1.67± 0.08 


13.63± 0.15 


1.55± 0.11 


21.31± 0.52 


2 


= 


1.36± 0.17 


18.05± 0.60 


1.63± 0.17 


11.75± 0.28 


1.61± 0.06 


11.75± 0.28 


1.57± 0.09 


19.28± 0.28 



growth in the number density of objects with temperature 
greater than about 0.5 keV. This corresponds to the later 
formation and propagation of the shocks in these models. 
In fact after z — 0.5 the situation becomes similar to that 
of the low-f2sM models, and the temperature distributions 
show little further evolution. 



4.5 Luminosity- Temperature relation 

In Figure 8 we present the distribution of emission- weighted 
temperature of the clusters as a function of the X-ray 
bolometric luminosity at three different redshifts: 2 = 1 
(crosses), z = 0.5 (open circles) and z — (filled cir- 
cles). For all the models there is a similar trend in the 
luminosity-temperature relation, even though the total num- 
ber of objects grows with the baryon fraction. We observe 
that a given temperature corresponds to higher luminosities 
in high-fisM models. This is related to the higher baryonic 
densities present in these cases. 

The luminosity-temperature distributions have been fit- 
ted by using a power-law relation of the form T = 10 b L x , 
where T is in keV and L x is in units of 10 40 erg s -1 . The re- 
sults, reported in Table 6, show that for all the models there 
is a little evolution with time, especially between z = 0.5 
and z — 0, where present-day observational data are avail- 
able. In particular the parameter r) is always inside the 1-a 
errorbar range, even if the models with smaller baryonic 
content tend to have a steeper slope. The normalization 6, 
instead, slightly decreases with increasing redshift and/or 
CIbm- This result is in qualitative agreement with the more 
recent observational analysis, which found no evolution of 



the temperature-luminosity distribution, at least for z < 0.5 
(Mushotzky & Scharf 1997). 

Our results can be directly compared to the observa- 
tional data. The hatched region in Figure 8 shows the L x -T 
relation (always with 1-a errorbars) obtained by Henry & 
Arnaud (1991) for clusters with luminosity larger than 10 44 
erg s~ . Even if the statistics are poor, the most luminous 
clusters for all the models are in good agreement with these 
data. The dotted line, instead, shows the fit from the com- 
bined sample of David et al. (1993), which contains clusters 
with lower luminosities (L x > 10 42 erg s _1 ). In this case the 
dispersion of the data around the fit (not explicitly reported 
in the original paper) is of the same order of that showed for 
the Henry & Arnaud (1991) results. The models with high 
baryonic fraction, even if they reproduce well the slope of the 
relation, have a lower normalization: at the same luminosity, 
the temperature is at least a factor 3 smaller than for the 
observations. The agreement is better for BM05 and BM10 
models which have a steeper (but still consistent) slope. 

Notice that it is not possible to present a comparison 
with the relation obtained by Mushotzky & Scharf (1997) 
because it refers to clusters with luminosities higher than 
those reached in our simulations. 



5 CONCLUSIONS 

In this work we have studied the evolution and the properties 
of X-ray clusters of galaxies in four different critical-density 
CDM-like models, in which the baryon fraction has been 
varied from Qbm = 0.05 to Qbm = 0.20. Models with a 
baryonic content larger than the predictions of the standard 
nucleosynthesis have been firstly considered by White et al. 
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Table 6. The parameters of the fits of the Luminosity- Temperature relation T = 10 b L^, computed at various redshifts for the different 
models. 



BM05 



BM10 



BM15 



BM20 







V 


b 


V 


b 


V 


b 


V 


b 


z 


= 1 


0.41±0.04 


-1.22±0.09 


0.34±0.03 


-1.31±0.07 


0.30±0.04 


-1.38±0.10 


0.31±0.04 


-1.52±0.11 


z 


= 0.7 


0.40±0.03 


-1.10±0.07 


0.33±0.03 


-1.20±0.09 


0.31±0.04 


-1.32±0.09 


0.30±0.03 


-1.42±0.08 


z 


= 0.5 


0.42±0.03 


-1.06±0.07 


0.35±0.03 


-1.17±0.09 


0.33±0.03 


-1.28±0.07 


0.30±0.03 


-1.35±0.07 


z 


= 0.2 


0.42±0.02 


-0.99±0.06 


0.36±0.02 


-1.11±0.06 


0.33±0.02 


-1.17±0.06 


0.31±0.03 


-1.25±0.08 


z 


= 


0.40±0.02 


-0.90±0.06 


0.37±0.02 


-1.05±0.05 


0.35±0.02 


-1.14±0.05 


0.32±0.02 


-1.20±0.06 



(1996) who found that they are in good agreement with a 
large set of observational data when coupled with a small 
tilt in the primordial spectrum (n w 0.8). Our results have 
proved to be useful in order to discriminate between the 
various models and to decide which of these models, if any, 
is compatible with observations. 

The mass function, the luminosity function, and the 
luminosity-temperature relation are the quantities that gave 
the most important hints on the properties of the models. 
The behaviour of these quantities is determined by the dy- 
namical evolution of the clusters. This is driven from the 
balance of two phenomena: the gravitational collapse of sin- 
gle objects and the merging of different structures. In the 
early stages of the evolution, the first effect tends to dom- 
inate, and the X-ray emission grows rapidly. The collapse 
of the baryonic matter is stopped by the formation of the 
shock. This rises strongly the pressure of the matter that 
finally is able to contrast the gravitational infall. Then the 
smaller virialized objects start to merge together, forming 
larger structures characterized by smoother density fields 
and hence by lower X-ray emission. 

The mass function presents the expected behaviour, 
with the cluster number density that, at z — decreases 
with increasing baryon fraction. This is mainly due to the 
different amount of power on small scales in the initial spec- 
tra of the density fluctuations. This result is obtained by 
identifying clusters only by their mass, without considering 
their luminosity. On the other hand, when clusters are se- 
lected by their total X-ray luminosity, the opposite trend 
is found: high-fisM models have the higher number density 
of X-ray clusters. This is essentially related to the presence 
of more baryons and so to the formation of higher baryonic 
density peaks that then leads to much higher X-luminosity, 
this quantity depending on the square of the baryonic den- 
sity itself. The X-ray emission depends also on the cluster 
temperatures but these are roughly the same for all the mod- 
els. 

Another difference is shown by the time evolution: while 
the number of clusters with large mass is a growing function 
of time for all the models, the abundance of luminous X- 
ray clusters starts to decrease at a some redshift, which is 
dependent on the cosmological model, being lower for high- 
!1m models. 

We compare the predictions of the four different cosmo- 
logical models to a series of observational results, mainly re- 
ferring to local (z « 0) datasets. By analysing the luminosity 
function and (more marginally) the luminosity-temperature 



relation we can conclude that the models with low bary- 
onic content (Qbm~ 0.05) are in better agreement with the 
data, while models with 15 or 20 per cent are well outside 
the 1-a errorbars. On the contrary the study of the clus- 
ter temperature distribution cannot distinguish between the 
various models. Our result is even more stringent if we ob- 
serve that our resolution tends to underestimate the lumi- 
nosities; consequently the differences of the high-fisM mod- 
els with observations would be even higher. Therefore the 
X-ray properties seem to exclude that the increase of the 
baryonic content can help to reconcile, in the framework of 
critical-density models, the cold dark matter scenario with 
the observations. 

In order to investigate if the differences found in the 
results for the various cosmological models are due to the 
different power spectrum normalization (i.e. as), we ran a 
second simulation for the BM15 model, but adopting in this 
case the same ag value of the BM05 model which we found 
previously to be the model in better agreement with the X- 
ray cluster data. As expected, in this new simulation with 
a higher normalization (hereafter BM15hn) the number Nm 
of clusters with mass greater than 1O 14 M0 is larger than in 
the BM15 model: we found N M = 28 and 72 at z = 0.5 
and z = 0, respectively. These results are comparable to 
those obtained for the BM05 model, confirming that the 
mass function depends mainly on the spectrum normaliza- 
tion. A similar comment can be done for the rms of the 
dark and baryonic matter and for the mean temperature 
computed on the cell scale. We also observe a larger for- 
mation of X-ray clusters in the BM15h„ model than in the 
BM15 model: at z = we found N 42 = 136 and N 43 = 38. 
The resulting luminosity functions, both bolometric and in 
the energy band [0.5-2] keV, are slightly higher than those 
of the BM15 model, increasing the difference with the ob- 
servational datasets. A stronger discrepancy with the data 
is also obtained when the luminosity-temperature relation 
is considered. All these results support the idea that an in- 
crease of the primordial baryonic fraction with respect to the 
standard values (approximately 5 per cent) produces X-ray 
clusters with properties in disagreement with the available 
data. Of course this is obtained by assuming a particular 
set of parameters, namely n = 0.8 and h — 0.5 and it would 
be interesting to know if the conclusion survives when these 
parameters are varied inside the range still allowing an ac- 
ceptable fit to other data (see White et al. 1996). 

Finally, another important feature which we found in 
our simulations is the very little redshift evolution of X- 
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ray cluster properties. In all of our models, the luminosity 
function, the temperature distribution and the luminosity- 
temperature relation for the simulated clusters are almost 
constant, particularly between z = 0.5 and z = 0, in good 
agreement with that which seems to emerge also from very 
recent observational data. 
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BM05 BM10 




Figure 1. a. The contour plots for the baryonic density Q BM in a slice of 64 X 64 X 0.5 h~ 3 Mpc 3 at z = for the four different models: 
BM05 (top left), BM10 (top right), BM15 (bottom left), BM20 (bottom right). The baryonic density is normalized to its mean density 
and the contour levels correspond to lO^ -3 '/ 4 , where i = 1,2, 15 
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Figure 2. * 

b. The same as Figure la but for the dark matter density Q DM which is normalized to its mean density. The contour levels correspond 

to lO^" 3 )/ 4 , where i = 1, 2, 14 
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log T 

Figure 3. * 

c. The same as Figure la but for the temperature T which is in units of Kelvin degrees. The contour levels correspond to 10 2 */ 3 , where 

j = l,2,...,12 
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Figure 4. * 

d. The same as Figure la but for the X-ray luminosity L x which is in units of 10 36 erg s _1 . The contour levels correspond to 10 28 / 3 , 

where i = 1, 2, 11 
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Figure 5. The contour plots for the baryonic density Q BM (first column), the dark matter density Q DM (second column), the temperature 
T (third column) and the X-ray luminosity L x (last column) in a slice of 20 X 20 X 0.5 h~ 3 Mpc 3 around the most luminous cluster in 
the BM05 simulation. The different rows show the redshift evolution: z = 1, z = 0.5, z = 0.2 and z = from the top to the bottom. The 
density of each component is normalized to its mean density while the temperature and the luminosity are in units of Kelvin degrees 
and 10 36 erg s _1 , respectively. The density contour levels correspond to lo' l ~ 3 '/ 4 while the temperature and the luminosity levels are 
10 2i / 3 , where i = 1,2, ... 
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Figure 6. The number of clusters Nm with mass > 1O 14 M0 as a function of the redshift z for the different models: BM05 (solid line), 
BM10 (dotted line), BM15 (short-dashed line) and BM20 (long-dashed line). 
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Figure 7. The number of clusters with luminosity L x > 10 42 (W42, upper panel) and L x > 10 43 erg s _1 (W43, bottom panel) 
as a function of the redshift 2 for the different models: BM05 (solid line), BM10 (dotted line), BM15 (short-dashed line) and BM20 
(long-dashed line). 
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Figure 8. The cluster bolometric luminosity function for the different models: BM05 (top left), BM10 (top right), BM15 (bottom left) 
and BM20 (bottom right). The different curves refer to various rcdshift: 2 = (solid line), z = 0.2 (dotted line), z = 0.5 (short-dashed 
line), z = 0.7 (long-dashed line), z = 1 (dottcd-dashed line). The hatched region shows the observational results (with l-cr errorbars) 
obtained by Ebeling et al. (1997). 
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Figure 9. The same as Figure 5 but in the [0.5-2] keV energy band. The vertically hatched region shows the observational results (with 
1-cr crrorbars) obtained by Ebeling ct al. (1997), while the horizontally hatched one refers to the Dc Grandi (1996) results. 
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Figure 10. The cluster temperature function for the different models: BM05 (top left), BM10 (top right), BM15 (bottom left) and 
BM20 (bottom right). The different curves refer to various rcdshift: 2 = (solid line), 2 = 0.2 (dotted line), z = 0.5 (short-dashed line), 
z = 0.7 (long-dashed line), 2 = 1 (dotted-dashcd line). The hatched region shows the observational results (with l-cr errorbars) obtained 
by Henry & Arnaud (1991). 



© 0000 RAS, MNRAS 000, 000-000 



20 Claudio Gheller, Ornella Pantano and Lauro Moscardini 



0.5 - 







o -0.5 



0.5 







n 1 1 r 



o • £ x * 



'oo o 

o 



t — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — rpp 



4 ^t^Vo 




x Z= 1 

o z=0.5 
• z = 



BM05 



i i i 



i i i i —a 



Q-- 



O 

• O 

*• -i • 





BM10 



i i i | i i i i U 



P.*' 



Q» jp TO 

••V^>o x x* 

O x 



« b o 



^x 



J I I I I I I I I I I I I L 



BM15 




BM20 



_3 I I I I I I I L_ 



-2 -1 

log L x (10 44 erg s" 1 ) 



-2 -1 

log L x (10 44 erg s" 1 ) 



Figure 11. Scatter plots of the bolometric X-ray luminosity L x and the emission-weighted temperature T of the clusters for the different 
models: BM05 (top left), BM10 (top right), BM15 (bottom left) and BM20 (bottom right). Different redshifts are displayed by different 
symbols: 2 = 1 (crosses), z = 0.5 (open circles) and z = (filled circles). The hatched region shows the observational results (with 
l-a errorbars) obtained by Henry & Arnaud (1991); the dotted line refers to the fit obtained from the combined sample of David et al. 
(1993). In this last case the errorbars are not explicitly reported, but they are of the same order of magnitude of the previous ones. 
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